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ABSTRACT 

We calculate the rapidity distribution and the total cross section of coherent and incoherent 
^ \ /9-production in the heavy ion ultraperipheral collisions (UPC) at ^/snn = 130 GeV using the gen- 

eralized vector dominance model(GVDM) and the Gribov-Glauber approach. We find the coherent 
cross section of p-production acoh = 490 mb compared to the acoh = 370 ± 170 ib 80 mb recently 
reported by the STAR collaboration at RHIC. The predicted cross section inside the acceptance of 
the experiment, \y\ < 1, agrees with the data within half a standard deviation. It is emphasized that 



. measurements of the rapidity distribution will provide a much more stringent test of the model. 

> : 
•1— I , 

^ ■ 1 Introduction 



Ultraperipheral collisions(UPC) of relativistic heavy ions at RHIC and LHC open a promising new 
avenue for experimental studies of the photon induced coherent and incoherent interactions with 
nuclei at high energies (see fl], I, I for the reviews and extensive lists of references). Really, the 
LHC heavy ion program will allow studies of photon-proton and photon - nucleus collisions at the 
energies exceeding by far those available now at HERA for 7 — p scattering. 

Hence, it is very important to check our basic understanding of the UPC processes using the re- 
actions which have smaller theoretical uncertainties on the level of the 7 A interactions. Recently the 
STAR collaboration released the first data on the cross section of the coherent p-meson production 
in gold-gold UPC at \/snn = 130 GeV This provides a first opportunity to check the basic 

features of the theoretical models and main approximations which include the Weizsacker- Williams 
(WW) approximation for the spectrum of the equivalent photons, an approximate procedure for 
removing collisions at small impact parameters where nuclei interact strongly, and the model for 
the vector meson production in the 'jA interactions. In the case of the p-meson production the basic 
process is understood much better than for other photoproduction processes. Hence, checking the 
theory for this case is especially important for proving that UPC could be used for learning new 
information about photon - nucleus interactions. 

Earlier we published predictions for the cross section of this process at higher energy ^/sly]\^ = 
200 GeV. Hence, a direct comparison of the STAR result with 0] is difficult. In this paper we perform 
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the analysis of the /9-production at ^/sNN = 130 GeV including the effects due to the cuts of the 
STAR experiment. 



2 Outline of the model 



Production of p-mesons in ultraperipheral heavy ion collisions can be expressed in the Weizsacker- 
Williams approximation H through the jA pA cross section 



dcr AiA2^AiA2V 
dy 



n\ {y)a^A-,^pA2 iv) + {-y)(j^A^^pA^ (-y)- 



(1) 



The quantity y = In^^ is the rapidity of the produced p meson and n'^{y) is the flux of photons 
with the energy uj^ = ^cQo emitted by one of nuclei (7c is the Lorentz factor for colliding nuclei, and 
qo is the photon momentum in the coordinate system of moving nucleus). The photoproduction 
cross section a^A^pAiv) can be calculated in the Glauber model |^] 



CTyA^pAiy) 



dt-r^\F^A-.pA{t)\ = -ra / t^^J 
p p i 



^ / dbe''^^''T{b) \ 

2TT J 
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Here = t± = tm, 
r(6) is the diffractive nuclear profile function 

T(b) = lim ^{h,z) 



-r-2- is longitudinal momentum transfer in 7 — p transition, and 



(3) 



To calculate the eikonal function $(6, z) the Glauber approach |1C] was combined with the gener 



alized vector dominance (GVD) model |11]. More properly such an approximation should be called 
the Gribov- Glauber modeljl^ because the space-time evolution of high energy processes is different 
in quantum mechanical models and in quantum field theory and therefore theoretical foundations 
for the high-energy model are different. In particular, in QCD in difference from quantum mechan- 
ics a high-energy projectile interacts with all nucleons at the same impact parameter almost at the 



same time [12|. Such formulae allow to extend the domain of applicability of the Glauber model to 
the description of high energy phenomena where inelastic (high multiplicity) particle production 
gives dominant contribution to the total cross section. As we are mostly interested in the accurate 
estimate of the coherent diffractive production of the vector meson p with Mp = 0.77 GeV, the 
spectrum of the higher corresponding hadronic states of M < 2 GeV can be approximated by 
one effective p'-meson with some reasonable fixed mass, say Mpi = \/3Mp||]. We want to draw 
attention that the value and sign of the p' contribution as taken from quenched GVDM is fitted 
to describe approximate Bjorken scaling in the case of highly virtual photon. Thus the model used 
in the paper correctly accounts for the fluctuations of cross section including color transparency 



phenomenon |14|. 

Then the GVD model comprises elementary amplitudes 

6 6 
f-yN-^pN = -ffpN-^pN + -7-fp'N^pN, 
Jp Jp' 
6 6 
ffN^p'N = -7-fp'N^p'N + -ffpN^p'N- 
Jp' Jp 



(4) 
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In the optical limit {A ^ 1), with accuracy 0{^/aem) the eikonal functions <I>p^p'(6, z) are determined 
by the solutions of the coupled two-channel equations 

2ikp-^<S>p{b, z) = U^A-.pA{b, z)e<C''' + 

^)^p(6, z) + UpA~.Mb, z)e<^"'^p,ib, z) (5) 

2ikp,—^p,{b, z) = UjA-*piA{b, z)e''\\ + 
+Up,A^p,Aib, z) + Up,A^pAib, z)e'<'^'''^p{b, z) (6) 

with the initial condition <I>p p/(6, — cxd) = 0. The exponential factors exp[iq^'' z] are responsible for 

the coherent length effect, i,j = 7, p, p', = ^ ' . The generalized Glauber -based optical 
potentials in the short-range approximation are given by the expression 

UiA^jA{b, z) = -4:7i-fiN^jN{0)g{b, z). (7) 

Here fiM^jNi^) are the forward elementary amplitudes, and Q{b, z) is the nuclear density normalized 
by the condition / dbdz gijj, z) = A. We calculated Q{b, z) in the Hartree-Fock-Skyrme (HFS) model 
which provided a very good(with an accuracy 2%) description of the global nuclear properties 
of spherical nuclei along the periodical table from carbon to uraniump^ and the shell momentum 



distributions in the high energy (p,2p)[16| and (e,e'p)[17| reactions. 

Following the simple suggestion of Ref . ||l^ which is quite reasonable in the case of light vector 
mesons, we fixed the elementary scattering amplitudes and coupling constants by relations 

fp'N^p'N = fpN-^pN, fpN^p'N = fp'N~^pN = —^fpN^pN, fp' = • fp, (8) 

with the value of the p-meson coupling constant /p/47r=2.01. The diagonal amplitude fpN^pN 



was taken in parameterization of the Landshoff-Donnachie model 1 18] while the value of parameter 
e = 0.18 was found from a best fit to the differential cross section of the p- meson photoproduction 
off lead at oo^ = 6.2 GeV and t± = 0.001 GeF^fl^. With all parameters fixed, we compared our 
calculations with all available data on p-meson photoproduction off nuclei at low and intermediate 



energies |19] and found a very good description of the absolute cross section and of the momentum 
transfer distributions (see Q). Hence, it was natural to expect that this model should provide a 
reliable parameter-free predictions for production of p-mesons in high energy heavy ion UPC. Note 
here that the inelastic shadowing effects which start to contribute at high energies still remain a 
few percent correction at energies < 100 GeV relevant for the STAR kinematics. For LHC energy 
range one should account for the blackening of interaction with nuclei. In this case cross section of 
inelastic diffraction in hadron- nucleus collisions should tend to 0. Also it leads to a suppression of 
the p' contribution to the cross section of the diffractive p meson photoproduction [^. 

3 Results and discussion 

The calculated momentum transfer distributions at the rapidity y = and the momentum transfer 
integrated rapidity distribution for gold-gold UPC at y/s^N = 130 GeV are presented in Figs.||a,b. 
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Figure 1: (a) Momentum transfer dependence of the coherent and incoherent p-meson production 
in AuAu UPC at \/snn = 130 GeV calculated in Generalized Glauber model (GGM). (b) Rapidity 
distributions for coherent p-meson production in the gold-gold UPC at ^/sattv = 130 GeV calculated 
in GGM. 
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Figure 2: Energy dependence of the cross sections for coherent /?- meson production in the gold-gold 
UPC calculated in the GGM and the STAR results: a) The dashed line is the total cross section in 
GGM, the solid line is the cross section calculated accounting for the STAR cut on the momentum 
transfer, star is the STAR cross section based on the Monte-Carlo extrapolation of the measured 
value to the full detector acceptance; b) Comparison of the GGM cross section in the interval of 
rapidities | y |< 1 with the value measured by the STAR. 
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Let us briefly comment on our estimate of the incoherent p-meson production cross section. The 
momentum transfer distribution (dashed hne in Fig.||a) is practicahy flat in the discussed t± range. 
The total incoherent cross section obtained by integration over the wide range of t± is dine = 120 
mb. To select the coherent production the cut t± < 0.02 GeV^ was used in the data analysis [^. 
Correspondingly, the calculated incoherent cross section for this region of t± is (Tj„c = 14 mb. Our 
calculations of incoherent production which are based on accounting for only the single elementary 
diffractive collision obviously present the lower limit. The residual nucleus will be weakly excited 
and can evaporate only one-two neutrons. The events A + A ^ p + xn + Ai + A2 were detected 
by the STAR and identifled as a two-stage process - coherent /O-production with the subsequent 
electromagnetic excitation and neutron decay of the colliding nuclei [^^. In particular, the cross 
section estimated by the STAR for the case when only one of the nuclei is excited and emits several 
neutrons is o"^^ = 95 it 60 it 25 mb. The momentum transfer distribution for these events is 
determined by the coherent production. Hence, it differs from that for incoherent events but in 
the region of very low t±_ it is hardly possible to separate them experimentally and obviously the 
measured cross section gn includes contribution of incoherent events on the level of 15%. 

The total rapidity-integrated cross section of coherent p-meson production calculated in the 
GVDM for the range of energies available at RHIC is shown in Fig.^ (dashed line). We find 
af^^ = 540 mb at ^snn = 130 GeV. The value al^l = 370 ± 170 ± 80 mb was obtained at 
this energy by the STAR from the data analysis at the low momentum transfer < 0.02 GeV^. 
Thus, before making a comparison we should take into account this cut. It leads to a reduction of 
the cross section by 10% (the solid line in Figj2|a). In our calculations we did not account for 
the t^-dependence of the elementary amplitudes which are rather flat in the considered range of 
energies and momentum transfers as compared to that for the nucleus form factor. So, in the region 
of integration important for our analysis neglect by this slope is a reasonable approximation but, 
nevertheless, an account of this effect would slightly reduce our estimate of the total cross section. 
Also we neglected a smearing due to the transverse momentum of photons and the interference of the 



production amplitudes from both nuclei [23|. This latter phenomenon results only in the narrow 
dip in the coherent t_|_-distribution at t± < 5 • 10~^ GeV'^. All these eff'ects does not influence 
noticeably the value of the tj_-integrated cross section but can be easy treated and taken into 
account in a more reflned analysis. Thus we flnd cr*^^ = 490 mb to be compared to the STAR value 
'^coh ~ "^"^^ -'- ^'''^ -'- Since our calculation does not have any free parameters, this can be 

considered as a reasonable agreement. 

At this point we would like to comment on the statement of @ that our prediction for coher- 
ent /9-production in gold-gold UPC at -y/sArAr=200 GeVp] is 50% higher than the value given by 



model in[21|. We already briefly explained in |g] that this discrepancy originates from a number of 



approximations made in the model of Ref.|21| which differs from the conventional Glauber model. 



The coherent photoproduction cross section was deflned in|21] by the expression: 



o,A^,A = '^''^^-"f = . j \Fm'dt. (9) 

—00 

where Fa (t) is the nuclear form factor (two-dimensional Fourier transformation of the parameterized 
nuclear density) and the forward photoproduction cross section was estimated using the vector 
dominance model and optical theorem 

'^'""-f = = ^-LipA). (10) 
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The total cross section of the pA interaction was found in 



1 21] using the formula: 



db 



1 — exp 



(^pN / Q{b,z)dz 



(11) 



Eq.(ll) presents the classical mechanics model with standard for this approach expression for the 
total cross section af^{pA). The quantum mechanics expression is given by the Gribov-Glauber 
model (here for simplicity we give the expression in the hmit of Re/Im=0) []: 







1 db 


1 — exp 



--apNTA{b) 



(12) 



In the Black Body Limit (uy at 
the classical mechanics {&- 



tot 



-> oo 



the total VA cross section estimated with the use of 

and the quantum mechanics (cr^^' = 2ttR\) expressions in 



Eq.(lO) differ by a factor of four. The difference for the case of the gold nucleus and reasonable value 
of the pN elementary cross section ~ 25 mb can be found using the simplified model of the 
nucleus of constant density qq 0.16 fm^'^ and radius Ra- With the value of radius Rau ~ 6.5 fm 
one can obtain the reasonable estimate of the ratio 



1 



2^2 „2d2 



'pN 



QqR^ 



1.55 



Hence, we find that due to use of the classical mechanics expression (^) instead of the Gribov- 
Glauber model expression (12) the total ^A — > pA cross section was underestimated in Ref.|21] 
by a factor ~ 2.5. Note, that in the range of the photon energies essential in photoproduction of 
p-mesons in UPC at RHIC the elementary pN cross section still weakly depends on the energy of 
/9-mesons. Hence, this factor weakly depends on ^/sNN■ 

Next important approximation made in Ref.||2l| in deriving formulae Eqs.(|9|) and (|lO| ) was the 

neglect by the coherence length effects. This requires q'u^^z = „ <^ 1. This neglect is not 

II -^IcQO 

justified because large longitudinal distances are essential in the diffractive p-photoproduction. The 
/9-meson can be formed far from the nucleus. Besides, the photon flux is large at small qo ^ Ra^^ i-e- 
in the region where the coherence length effect is important. We estimated that the cross section for 
AuAu pAuAu at the energy ^stvat = 200 GeV is overestimated by a factor « 1.5 if one neglects 
the coherence length effect. The coherence length effect becomes more essential with decrease of 
the energy. As a result at ^/sNN = 130 GeV this effect suppresses the cross section by a factor 2. 
On the contrary, at much higher energies where the coherence length is very large (for example, at 
LHC) this effect will be small. 

Note in passing that the calculations in Ref. |pT|| were performed neglecting the real part of the 
elementary pN amplitude which we accounted for using the Landshoff-Donnachie parameterization. 
In the high energy domain, for example, in the region of the central rapidities at RHIC the real part 
of the pN amplitude is negligible but one should account for RefpN^pN at the edges of rapidity 
distribution which correspond to the photoproduction at intermediate energies of photons. Since 
the contribution of this region to the total cross section is enhanced by the high photon flux, the 

^In the Appendix we demonstrate how the model used in Ref.j2l| but with the correct high energy expression for 
the total cross section can be obtained from the Gribov-Glauber based Generalized Vector Dominance Model and 
what essential approximations have to be done on the way. 
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total cross section of the coherent p production at ^JsNN = 130 GeV would be underestimated by 
« 10%. 

Thus we emphasize that a number of approximations were made in the model of the Ref. pl[| 
which are inconsistent with the conventional Glauber model and their interplay depends on the 
photon energy. Thus the estimates obtained within this model contain systematic uncertainties. 
We want to draw attention that the model we have used is well theoretically justified. It cor- 
rectly calculates the nuclear form factor in the coherent photoproduction. We checked |^] that this 
model provided a very good description of the coherent /o-photoproduction off nuclei at low and 
intermediate energies along the periodical table without any free parameters. 

At present the comparison of our predictions with STAR data is still preliminary because the 
experimental errors are too large and there are few points in the procedure of the data analysis 
which should be discussed. Really the acceptance of STAR is very strongly y dependent being 
maximal at y = and going to zero at | y |= 1 (Fig. 3a in [^) while the theoretical distribution is 
expected to have a double bump shape, see Fig.|l]b, which is simply due to the symmetry of collision 
and the interplay of the energy dependence of the photon flux and pA cross section. 

Due to the acceptance conditions AuAu pAuAu events were detected in the range of rapidities 
\y\ ^ 1) while the cross section reported in Q is corrected for the | y |> 1 using the Monte-Carlo 
extrapolation based on the model in which ^ 

R = aiTr/cr\y\<i = 2.7. (13) 

Thus the cross section in the region | y |< 1 is 

a|^l<i = 140 ± 60 ± 30 mb. (14) 



The errors for this cut are scaled accordingly. However the model |21] differs from our in many 
aspects, described above. In particular, the neglect by the coherence length effect leads to a signif- 
icant modification of the value of factor R. The rapidity interval | y |< 1 corresponds high energies 
of the photon where this effect is negligible. On the other hand due to the neglect by the coherent 
length effects at the edges of the rapidity distribution the relative contribution of photoproduction 
cross section at | y | > 1 was overestimated in Ref . [^l| . Hence their value of R is lower than the one 
we obtain from our rapidity distribution which gives R = 3.2, and the cross section (^\y\<i = 170 
mb. It is this value of the cross section which should be compared (Fig.|2|b) to the experimentally 



measured cross section (14). Obviously a more detailed comparison would require a detailed study 



of the sensitivity of the analysis to the assumed y-distribution. 



4 Conclusions 

We demonstrate that the cross section of coherent p-meson production in high energy heavy ion 
UPC calculated within the GVDM is in a good agreement with experimental results of the STAR 
collaboration. More stringent test will involve comparison of the model predictions with the cross 
section and the rapidity distributions of the p-meson production measured with higher precision. 

We thank A. Baltz, B. Kopeliovich, P. Yepes, and especially S. Klein for very useful discussions. 
This study was supported in part by GIF, DOE and CRDF 

^We thank S.Klein for the discussion of this issue. 
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5 Appendix 



For the completeness let us discuss how classical mechanical formulae used in Ref.|21] arise from the 
quantum mechanical ones. Since the Vector Dominance model was used in [^] we should neglect 
by the nondiagonal p ^ p' transitions in Eq.(^. Thus we need to put e = in Eq.(^ instead 
of e = 0.18 which is used in our calculations. If we keep all other parameters of the model fixed, 
we pay for such a reduction by the increase of the p-photoproduction cross section by a factor 
« (1 + 2e/V3) 1.2. With e = equations (|) and (|) become decoupled and the solution of 
Eq. (^) gives the eikonal function of the p-meson 



^p{b,z) 



1 2h-p ! UpN^pN{b,z')dz 



2z/l/i 



dz'U./N-*pNib,z')e'^' 



(15) 



Using Eq.(|2|) and Eq.(|^) and standard expression for the elementary amplitude 



fpNiO) 



ika 



'pN 



1 — ifipN] PpN 



RefpNiO) 



47T - /m/p^(0)' 
we obtain the amplitude of the p-meson photoproduction off the nucleus in the optical limit of the 



standard Glauber plus the Vector Dominance model 1 22 1 



F^A^pA = fjN^pNiO) J dhe'^^ 




dz' Q{b, z')e''\\ 



, -^[l-i/3pjv] / e{b,z")dz" 



(16) 



Note that accounting for the real part of the pN amplitude leads to appearance of the phase 

00 _ 

factor exp[i(5pN J Q{b, z')dz'] which is similar to that describing the coherence length effect and 



which is important in the same energy domain. Following the assumptions of Ref.|21] where both 
the coherence length effect and real part of the pN amplitude were neglected we remove these 
exponential factor from Eq.(|l6|). Then in the limit of the purely imaginary elementary pN amplitude 
we obtain 



F^A^pA = f-yN~>pNm J dbe''^^-' J dz'gib,z')e 
-00 



/ e(b,z")dz" 



f^N^pN{0) 



O'pN 

_ f-/N^pNiO) 
(TpN 



■2- J dbe^'i^-^ 




dz —exp 
dz 



■2- J dbe'^^-^ 




1 — exp 



7pN 



CpN 



g{b, z)dz' 



Q{b, z)dz 



(17) 



(18) 



Now using the vector dominance relation f^]\[^p]\f = -fzfpN^pN we can write the formula for the 

J p 

forward — > pA cross section in the optical limit of the Glauber +VD model 



da^A^pA{t = 0) _ vr 



dt 



— \F. 



■yA^pA 



it 



a 



em 2 



^Pp 



(19) 
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where 



apA = 2- db 



1 — exp 



Q{b, z)dz 



(20) 



Thus in the Glauber+VD model we have got the expression for the forward cross section of the 
photoproduction coinciding with Eq.(p!0|) used in calculations performed in Ref.|21] but with high 
energy quantum mechanics formula for the total cross section. 
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